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ABSTRACT 

We study the Kennicutt-Schmidt relation between average star formation rate and average cold 
gas surface density in the HI dominated ISM of nearby spiral and dwarf irregular galaxies. We 
divide galaxies into grid cells varying from sub-kpc to tens of kpc in size. Grid-cell measure¬ 
ments of low SFRs using Ha emission can be biased and scatter may be introduced because of 
non-uniform sampling of the IMF or because of stochastically varying star formation. In order 
to alleviate these issues, we use far-ultraviolet emission to trace SFR, and we sum up the fluxes 
from different bins with the same gas surface density to calculate the average Ssfr at a given 
value of Tigas- We study the resulting Kennicutt-Schmidt relation in 400 pc, 1 kpc and 10 kpc 
scale grids in nearby massive spirals and in 400 pc scale grids in nearby faint dwarf irregulars. 
We And a relation with a power law slope of 1.5 in the Hl-dominated regions for both kinds 
of galaxies. The relation is offset towards longer gas consumption timescales compared to the 
molecular hydrogen dominated centres of spirals, but the offset is an order-of-magnitude less 
than that quoted by earlier studies. Our results lead to the surprising conclusion that conver¬ 
sion of gas to stars is independent of metallicity in the HI dominated regions of star-forming 
galaxies. Our observed relations are better fit by a model of star formation based on thermal 
and hydrostatic equilibrium in the ISM, in which feedback driven turbulence sets the thermal 
pressure. 
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1 INTRODUCTION 

How gas is converted into stars is an important question crucial 
to our understanding of galaxy formation an d evoluti o n. An em¬ 
pirical way of quantifying th is was given by ISchmidtl ( Il959h and 
extended by iKennicutn d 1 998h . The Kennicutt-Schmidt (K-S) rela¬ 
tion is a power law relation between Esfr and Egas averaged over 
the star-forming disks of nearby galaxies. Insight into the physical 
processes regulating star formation has been sought by trying to 
model the K-S relation J&urnholz. McKee & Tumlinso'iil f2009bl : 

lOstriker. McKee & LerovIbolOh using physical principles as well 

as empirical clues. 

Only recently have surveys of nearby galaxies reached suffi¬ 
cient sensitivities at all requisite wavelengths to study the K-S re¬ 
lation on sub-k pc scales. Example s include the study o f nearby spi¬ 
ral gala xies bv lBigiel et al.l ( l2008h . In a similar study. Ilxrov et al.l 
( l2013ah showed that a linear K-S relation exists between the sur¬ 
face density of molecular gas and Esfr- This is not unexpected 
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as the molecular phase of hydrogen traces the molecular clouds 
where star formation actually occurs. A linear K-S relation be¬ 
tween Ehj and Esfr has also been shown to hold in the atomic 
hydrogen dominated out er discs of nearby spirals iSchruba et al.l 
l201lh]Bigieletal.l l l20ldl showed that in the very outskirts of spiral 
discs (outside R 25 ), where molecular hydrogen cannot be detected 
and atomic hydrogen (HI) dominates, there is a clear relationship 
between FUV and HI column density, with the star formation ef¬ 
ficiency increasing with increasing HI column. The inferred gas 
consumption timescales in these regions were well above a Hub¬ 
ble time (~ 10^^ years). Similar results were found in recent re- 
so lved sub-kpc studies of HI and sta r formation in dwarf galaxies 
bv iRovchowdhurv et alj ( l2009l . [20llh . 

One crucial factor to take into account when determining the 
relation between gas and star formation in the HI dominated regime 
is the low star formation rate surface densities in this regime. 
Very deep imaging data is required to measure fluxes accurately. 
In addition, since direct tracers of star formation primarily trace 
young high mass stars (> 16 Mq for Ha, > 6 Mq for FUV), the 
high mass end of the initial mass function (IMF) will be sampled 
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stochastically for low SFRs resulting in considerable fluctuation in 
the luminosity of the tracer across an ensemble of systems with the 
same SFR. Bursty SF histories will also result in errors when us¬ 
ing calibrations whi ch assume continuous star formation duri ng the 
previous Myr or so. Ida Silva. Fumagalli & Krumholzl bOldl) have 
simulated the effects of stochasticity and bursty SF history on the 
measurement of SFRs using different tracers. They show that that 
a bursty SF history leads to the peak of the distribution of SFRs 
measured for an ensemble of simulated systems to be offset on the 
lower side compared to the true SFR put in as input, and this offset 
increases with decreasing SFR. Stochasticity in the sampling of the 
IMF leads to the scatter in measured SFRs increasing wi th decreas¬ 
ing SFR. In a recent study, iRovchowdhurv et"^ l l20I4h computed 
a relation between SFR and atomic gas surface density in nearby 
dwarf irregulars averaged over their star-forming discs, and found 
that the Egas,atomic based K-S relation is not steep, but almost lin¬ 
ear. The implied efficiency of star formation was low in HI dom¬ 
inated dwarf galaxies, but it did not decrease as a function of HI 
surface density as found bv iBigiel et alj l l2010tl . A more systematic 
study of the relation between Esfr and Egas in the HI dominated 
regime is therefore now in order. Here we study the relation be¬ 
tween Esfr and Egas in nearby T* spirals and dwarf galaxies 
over different scales, carefully accounting for how we measure low 
SFRs. We start with kpc sized regions and zoom in further to study 
trends such as changes in scatter in the Esfr - Egas plane. We 
also study how the K-S relation changes for few tens of kpc sized 
regions in a sample of Lt, spirals with very large Hl-to-optical 
size ratios. 

2 THE SAMPLES 
2.1 THINGS spirals 

The first sample is of spiral galaxies from The HI Nearby Galaxy 
Survey (THINGS. IWalter et al.ll200il . THINGS is a high spatial 
and velocity resolution interferometric survey of HI in 34 nearby 
galaxies using the NRAO Very Large Array (VLA). The sample 
spans a wide range in luminosity and mass, and included a hand¬ 
ful of dwarf galaxies. The THINGS galaxies used in the present 
study have CO detect ions from the HE RA CO Line Extragalactic 
Survey (HERACLES iLerov et^l2009h . Spitzer 24 data, and 
low inc lination angles. This is the same sample used bv iLerov et al.l 
( l2013tJ) to derive the linear K-S relation between Esfr and Eh 2 • 
Our sa mple is also part o f the study of the very outskirts of spi¬ 
rals by dSigiel et aklboioh . All galaxies have GALEX FUV obser¬ 
vations which in combination with Spitzer 24 p,m data is used to 
estimate the i r SERs . We exclude one galaxy (NGC 3184) from the 
iLerov et^ ( l2013all sample whose GALEX observations are from 
the low exposure time All-sky Imaging Survey (AIS) rather than 
being from the Medium Imaging Survey (MIS) or Guest Investiga¬ 
tor data. Another galaxy (NGC 2841) was excluded for not having 
CO data sensitive enough to trace the full disk, especially the cen¬ 
tral regions where there is high star formation but a prominent hole 
in the HI distribution. Tracing the CO in the central regions of our 
sample galaxies is important because as mentioned in Section [T2l 
we use CO emission to exclude regions where HI does not domi¬ 
nate the gas in the ISM. 

The final sample of THINGS spiral galaxies are listed in Ta¬ 
ble □ Columns (1) gives the galaxy name, (2) the mass in stars 
(from ILerov et al.ir2013all . (3) the mass in HI (fr om I Walter et al. 
l2008h . (4) the distance to the galaxy taken from iKennicutt et^ 


Table 1. THINGS spiral galaxy sample 


Galaxy 

M* 

(10“ Mo) 

Mhi 

(10® Mo) 

D 

(Mpc) 

CO res. 
(kpc) 

HI cube 
400 pc 

NGC 0628 

1.0 

38.0 

7.2 

0.46 

NA 

NGC 2903 

2.5 

43.5 

8.9 

0.57 


NGC 3198 

1.0 

101.7 

14.1 

0.91 


NGC 3351 

1.3 

11.9 

9.3 

0.60 

NA 

NGC 3521 

5.0 

80.2 

11.2 

0.72 

RO 

NGC 3627 

3.2 

8.18 

9.4 

0.61 

NA 

NGC 4736 

1.6 

4.00 

4.7 

0.30 


NGC 5055 

3.2 

91.0 

7.9 

0.51 


NGC 5194 

3.2 

25.4 

7.9 

0.52 

NA 

NGC 5457 

2.5 

141.7 

6.7 

0.43 


NGC 7331 

6.3 

91.3 

14.5 

0.94 

NA 


( I 2 OIII) for all galaxies except NGC 5194 dWalter et alj|2008l) . Col¬ 
umn (5) gives the resolution in kpc corresponding to the CO res¬ 
olution 01-^13 arcsec. We require HI maps at linear resolutions 
of ~400 pc (see Section O, and seven of the galaxies had maps 
with average beam sizes comparable to 400 pc depending on the 
weighting scheme chosen: natural (NA) or robust (RO), and these 
are listed in column (6). 


2.2 EIGGS dwarf irregulars 

Our sample of dwarf irregular galaxies is derived from the Fain t 
Irregular galaxy GMRT Survey (FIGGS, iBegum et all l2008hll . 
FIGGS is a high spatial and velocity resolution interferometric sur¬ 
vey of HI in extremely faint (Mb > —14.5) gas rich (single dish 
HI flux > 1 Jy km s“^) dwarf galaxies within a volume of ~11 
Mpc. From the FIGGS sample we choose the galaxies which have 
archival GALEX far-ultraviolet (FUV) data, which can be used for 
estimating their SFRs. Once again galaxies observed as part of the 
low sensitivity AIS are excluded. We also exclude galaxies with 
HI distributions that are morphologically of fset from their stellar 
disc as found bv iRovchowdhurv et alj ( l2014l) . The FIGGS galaxies 
are highly neutral hydrogen dominated and have a typical gas-to- 
total baryon ratio of 0.7, and h ave extended HI disc s with the ra¬ 
tio of Hl-to-optical diameter ~2 l lBegum et al.l20085) . The FIGGS 
galaxies have no detectable CO emission which can be used to trace 
molecular hydrogen. 

The properties of the FIGGS galaxies used in this study are 
given in Table Column (1) gives the galaxy name. Some of the 
sample galaxies have available archival data from Spitzer based 
surveys of local galaxies. These galaxies are indicated with a ‘S’ 
in superscript after their names, followed by the reference to the 
original survey pap ers. Column (2) gives the mass in stars (see 
iBegum et al.ll2008al. for how the values were calculated), (3) the 
mass in HI (f rom IBegum et al] l2008hl) . (4) the distances to the 
galaxies (from lK^achentsev et al.ll2013l) . (5) the metallicity of the 
galaxy. For only a handful of the sample galaxies (indicated with 
the reference number as superscript) abundance measurements are 
available, fr om which the metallicities g iven in column (5) are cal¬ 
culated (see IRovchowdhurv et al.ll2014L for details). For the rest of 
the galaxies the metallicity tabulated is an estim ate using the lumi¬ 
nosity (Mb) - metallicity relation for dis from lEkta & Chengaluil 
( l2010h . These metallicity values are used to provide a correction 
factor to the measured SFR with the method described in Sec¬ 
tion [53] 
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Table 2. FIGGS dwarf galaxy sample 


Galaxy 

M* 

(10® Mq) 

Mhi 

(10® Mo) 

D 

(Mpc) 

Z 

(Zq) 

And IV 


205.19 

6.3 

0.06 

DDO 226 ®’1 

35.5 

25.95 

4.9 

0.12 

UGC 685®’! 

61.3 

56.15 

4.5 

0.20® 

KK 65 

74.3 

34.38 

7.62 

0.12 

UGC4115 

42.3 

285.53 

7.5 

0.12 

UGC 4459®’! 

20.1 

64.2 

3.56 

0.13® 

UGC5456®’! 

64.6 

58.96 

5.6 

0.16 

UGC 6456 ®’2 

43.0 

43.89 

4.3 

0.10“ 

NGC3741®’! 

14.0 

130.0 

3.0 

0.09® 

UGC 7242®’! 

31.2 

45.75 

5.4 

0.11 

CGCG 269-049®’! 

10.9 

26.4 

4.9 

0.05® 

UGC 7298 

10.8 

21.6 

4.21 

0.06 

DDO 125®’! 

69.6 

31.87 

2.5 

0.12 

UGC 7605®’! 

27.7 

22.29 

4.43 

0.10 

UGC 8215 

5.5 

21.41 

4.5 

0.06 

DDO 167 

11.5 

14.51 

4.2 

0.07 

UGC 8638®’! 

38.4 

13.76 

4.27 

0.10 

DDO 181®’! 

15.3 

27.55 

3.1 

0.14® 

DDO 183®’! 

9.6 

25.90 

3.24 

0.09 

UGC 8833 

7.8 

15.16 

3.2 

0.07 

DDO 187®’! 

10.1 

16.30 

2.5 

0.11® 

KKH98 

2.8 

6.46 

2.5 

0.04 


References- TjPale et a^j2009|h_2jEngelbracht etj^l f20o3) : 3: 

iMarble et alJ ilOldli 4lMou^kas & KemiicutiniOO^. 


2.3 Bluedisk spirals 

The final sample used in this study comes from a survey of 
11 > log (M^/Mq) > 10, 0.01 < z < 0.03 galaxies with un¬ 
usually high HI mass fractions - 0.6 dex hig her than the median re - 
lation between M(HI)/M* and M* found bv ICatinella et al.l | |201(]|) . 
Since these galaxies were chosen based on the bluer g — i colour 
in the outer optic al disc compared to the inner optical disc (see 
IWang et^l2013L for the details of the survey), they are termed 
‘Bluedisks’. These galaxies have HI discs extending much beyond 
their optical discs. An equivalently large sample of spirals with 
stellar masses and mass surface densities comparable to the pri¬ 
mary ‘Bluedisk’ sample were also observed in HI, and these form 
the ‘control’ galaxy sample. We exclude the galaxies identified as 
multi-source systems from our sample. 

The ‘Bluedisk’ and ’control’ galaxies do not yet have CO mea¬ 
surements, but here we only focus on the Hl-dominated regime by 
excluding the central regions (within R25) of the galaxies where 
Egas is likely to be dominated by molecular hydrogen. GALEX 
FUV data combined with WISE 22 gm data is used as the main 
tracer of SFR. For 27 of 42 sample galaxies, the source of the FUV 
data is AIS. We retain these galaxies for the study, because the re¬ 
gions over which we are calculating Esfr for this sample corre¬ 
spond to ~100 kpc“^ (~20 arcsec beams), so the signal-to-noise 
for the total FUV flux in each such region is large enough even for 
AIS data to be usable. 

Some properties of the ‘Bluedisk’ sample are listed in Table[3 
The upper section of the table lists the ’Bluedisk’ main sample, 
while the lower section of the table lists the ‘control’ sample. Col¬ 
umn (1) gives the galaxy identifier, (2) the derived stellar mass, 
(3) the derived HI mass, (4) the distance to the galaxy in Mpc, (5) 
the measured axial ratio in the r band (used when calculating the 
surface densities as described in Section [33] and given here due to 


Table 3. Bluedisk spiral galaxy sample 


ID 

M* 

Mhi 

D 

b/a 

HI beam 


(lo!® M0) 

(10® Mq) 

(Mpc) 


(kpc X kpc) 

Main 






1 

2.69 

12.30 

124.23 

0.64 

16.2x10.2 

2 

4.17 

8.71 

110.49 

0.41 

13.9X 9.2 

3 

2.63 

7.76 

128.88 

0.86 

16.1x10.8 

4 

3.72 

18.20 

115.06 

0.39 

17.5X 9.1 

5 

2.09 

15.85 

110.56 

0.89 

14.6X 9.4 

6 

6.92 

20.42 

110.55 

0.52 

15.lx 9.1 

8 

1.91 

15.85 

119.64 

0.74 

13.8x10.3 

12 

4.27 

13.18 

115.05 

0.45 

10.6X 9.1 

14 

6.17 

13.80 

105.97 

0.49 

11.3X 9.0 

15 

6.61 

24.55 

96.85 

0.73 

13.7X 7.9 

16 

2.04 

10.00 

119.70 

0.68 

12.6x10.2 

17 

4.90 

21.88 

124.21 

0.68 

12.8x11.1 

18 

2.19 

11.75 

101.46 

0.45 

12.5X 9.8 

19 

1.78 

14.79 

115.06 

0.81 

14.2X 8.9 

20 

1.62 

16.22 

115.12 

0.53 

14.lx 9.0 

21 

2.57 

7.59 

133.45 

0.60 

17.3x10.7 

22 

6.61 

10.00 

124.31 

0.28 

16.0X 9.6 

24 

5.13 

16.22 

133.42 

0.48 

13.0x11.4 

26 

2.04 

3.16 

101.45 

0.39 

12.4X 8.5 

30 

2.45 

11.75 

124.27 

0.26 

16.0x10.0 

35 

4.07 

8.51 

106.01 

0.36 

12.3X 9.1 

47 

3.39 

10.23 

105.99 

0.60 

15.3X 7.9 

50 

7.41 

6.31 

124.24 

0.29 

11.8x10.2 

Control 






9 

5.75 

8.13 

124.21 

0.45 

13.5X 9.1 

10 

8.13 

10.47 

128.91 

0.89 

12.8x11.1 

11 

4.37 

5.13 

105.99 

0.61 

11.2X 8.4 

23 

5.62 

8.71 

133.41 

0.43 

15.2X 9.7 

25 

3.47 

5.75 

124.21 

0.78 

12.0x10.2 

27 

2.00 

1.86 

128.84 

0.75 

17.8x10.0 

28 

3.47 

1.32 

83.33 

0.86 

9.0x 6.9 

32 

1.95 

3.63 

119.67 

0.44 

17.0X 9.4 

33 

5.50 

1.58 

119.67 

0.57 

13.1x10.4 

36 

2.19 

2.45 

110.53 

0.65 

11.3X 9.4 

37 

2.45 

5.01 

115.10 

0.69 

11.5X 9.8 

38 

5.62 

0.56 

110.48 

0.73 

10.9X 9.6 

39 

6.92 

1.38 

105.98 

0.75 

13.9X 8.3 

40 

2.29 

4.47 

119.70 

0.84 

17.1 X 9.9 

42 

6.31 

3.72 

124.28 

0.57 

13.6x10.7 

43 

1.82 

3.89 

119.67 

0.73 

15.0X 9.9 

44 

6.03 

1.91 

124.25 

0.64 

14.5 X 9.9 

45 

4.79 

4.07 

115.09 

0.39 

11.9X 9.6 

49 

2.63 

5.50 

133.42 

0.62 

20.6x10.9 

changes compared to values tabulatet 

inIWang 

et al 

boLlh. and 161 


the siz e of the HI beam in kpc. All values are taken from Wang et ^ 
( l2013h . 


3 ANALYSIS 
3.1 Overview 

We derive relations between Esfr and Egas or Egas,atomic in the 
Hl-dominated regime for the galaxies in our three samples. Here 
Elgas,atomic refers to the atomic gas surface density corrected for 
the presence of helium by multiplying the HI surface density by 
a factor of 1.34, and Egas refers to the total gas surface density 
including the contribution from H 2 . 
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We begin with (i) THINGS spiral galaxies at the resolution of 
the CO maps (mean and median resolution of ~600 pc) by aver¬ 
aging over 1 kpc sized regions. Next, we decrease the size of the 
regions over which the surface densities are compared. We consider 
(ii) THINGS spiral galaxies at a resolution of ~400 pc and, ( iii) 
FIGGS dwarf irregular galaxies at a resolution of ~400 pc. 400 pc 
is a natural and convenient choice of resolution for the T, spiral 
sample, because 7 of the 12 THINGS galaxies are observed in HI 
at a comparable native resolution. The resolution of the Spitzer 24 
p,m data is ~6", i.e. comparable to the HI. In the THINGS spiral 
sample, we choose Hl-dominated regions by referring to the CO 
maps and checking whether there is any detectable CO emission. 
This procedure ensures that not only do we choose choose regions 
from the outer disks of the galaxies, but also from the inter-arm 
regions of these spiral galaxies. We consider all regions in dwarfs, 
including the very central regions, as being dominated by HI. 

Finally we study the relation between Esfr and Egas at large 
(~10 kpc) scales using the ‘Bluedisk’ main and control spiral 
galaxies at the resolution of their HI maps. In this sample, all re¬ 
gions whose centres are within a distance of R25 from the optical 
centre of the galaxy are not considered. Once again, this is done in 
order to restrict the analysis to regions where the HI gas is likely 
to dominate. For all datasets, fluxes are extracted from matched 
square grid cells in the gas maps and the UV/IR images. Adjacent 
regions overlap by half the length of each side in order to ensure 
Nyquist sampling. 


3.2 Measuring surface densities of gas 

For the THINGS / HERACLES galaxies the flux in CO maps is 
converted to Eh 2 using the Galactic CO-to-H 2 conversion factor 
of 4.35 Mq pc“^ (K km The limiting E gas for HERA¬ 

CLES CO data is ~3 M© pc“^ l lLerov et alj 20091) . Only Egas val¬ 
ues above this limit are considered. 

HI maps are smoothed to the resolution of CO maps for the 
1 kpc scale study. For all other datasets HI maps are used at their 
native resolutions. The flux density in HI is converted to column 
density using the standard transformation for emission from an op¬ 
tically thin medium. T he limiting Egas.ato mic for THINGS galax¬ 
ies is ~0.5 Mq pc“^ dWalter et alj 200 Sh . and only regions with 
Egas,atomic greater than this value are considered. The 3cr HI col¬ 
umn density sensitivity threshold for FIGGS galaxies varies from 
0.4to5MQ pc“^ with a median of 1.4 Mq pc“^, and only regions 
in FIGGS galaxies with Egas,atomic greater than the corresponding 
threshold are used in our study. The 3a HI column density sensitiv¬ 
ity threshold s correspond i ng to the 2cr values given in Column (7) 
of Table 2 in IWang et alj ( 1201 3h are used as cut-off for ‘Bluedisk’ 
sample galaxies and only regions with higher Egas,atomic are con¬ 
sidered. 


3.3 Measuring surface densities of SFR 

Ida Silva. Fumagalli & Krumholzl | |2014|) have shown that Hq be¬ 
comes a highly unreliable tracer of SFR at low SFR. Therefore, we 
use FUV emission from GALEX (resolution ~4.5”). and use mid- 
infrared 24 pLxa emission from Spitzer (resolution ~6") or 22 p.m 
emission from WISE (resolution ~12") to correct our SFRs for at¬ 
tenuation due to dust. Foreground stars and background galaxies 
are identified and masked in the FUV and 24 pra maps of THINGS 
and FIGGS galaxies. For the THIN GS galaxies, this is done by 
applying the colour based masks of iMunoz-Mateos et alj l l2009h 


followed by visual inspection. For the FIGGS galaxies, the mask¬ 
ing is done by visual inspection. The masking of the WISE 22 pm 
images of the ‘Blued i sk’ sa mple galaxies is done as described in 
iHuang & Kauffmantil j2014h . The PSF of Spitzer 24 pm maps, and 
to a lesser extent that of the GALEA F UV maps have extend ed side- 
lobes. We use the kernels provided in lAniano et alj | |201 ih to con¬ 
vert to Gaussian beams close to the final required beam (that of the 
CO or HI map) for the 24 pm and FUV maps, and in the next step 
we smooth the Gaussian beam to match that of the CO or HI map. 
The WISE 22 pm maps are directly smoothed to the resolution of 
the corresponding HI map. Fluxes are then extracted from regions 
of requisite size in these smoothed maps. 

The conver sion between FUV luminosity and SFR is given by 

iHao et al I ItoIII) : 


log 


SFR 
Mq yr-1 


log 


Hz ergs 8“^ Hz“i 


43.35 (1) 


We account for the energy from star formation re-radiat ed at infra¬ 
red w avelengths due to dust using mid-infrared data iHao et all 

ll201lh : 


LpUV.corr = LpuV.obs + 3.89 L 25 /im 


( 2 ) 


which is fed into equation Q] to obtain the SFR. The calibration is 
valid for either Salpeter or Kroupa IMFs with stellar mases in the 
range of 0.1-100 Mq and solar metallicity. Only fluxes 3 times 
higher than the corresponding noise levels for Spitzer 24 pm and 
WISE 22 pm data are considered, as they are only additive correc¬ 
tions to the FUV flux. The a bove calibration i s also applicable to 
dwarf galaxies as discussed in lPe Looze et al.l ( l2014h . 

For FIGGS sample galaxies with 24 pm detections, we com¬ 
pare the detectable infrared fluxes in 400 pc regions to the FUV lu¬ 
minosity in the same regions. A bivariate regression fit to the data is 
then used to estimate the 24 pm flux in a region of a galaxy without 
Spitzer 24 pm measurements: 


log 


F24(im 

MJy st“i 


1.67 log 


Lpuv 

ergs Hz“i 


39.8 (3) 


This estimate is only derived for regions with sufficiently high FUV 
luminosity (> 24 ergs s“^ Hz“^). 

To estimate the error on the measured SFRs we add the fol¬ 
lowing terms in quadrature: measurement error for FUV and 24 
pm fluxes (as determined from the ‘relative response’ and ‘weight’ 
maps respectively), 10% flux calibration error for GALEX FUV 
data, 5% flux calibration error for Spitzer 24 pm data, and 50% 
error to account for the uncertainty in the SFR calibration caused 
by variations in t he IM F and star formation history following 
lLerovetalJ ( l20lll2013ah . 

The measured FUV flux used in equation also requires a 
correction if star formation is happening at very low metallicities 
as is the case for the FIG GS sample dwarfs (viz. Z ^ 0.1 Z q as 
can be seen from Tablel^. lRaiter. Schaerer & FosburvI (1201 Ol) used 
evolutionary synthesis models using a Salpeter IMF and constant 
star formation for the last 10® years to estimate emergent fluxes 
in sub-solar metallicity environments and found that FUV fluxes 
increase by ~ 11%, 19%, 27%, 32% for metallicities of 0.4, 0.2, 
0.05, 0.02 times solar. For ea ch FIGGS sample galaxy we do a 
linear interpolation between the lRaiter. Schaerer & FosburvI ( 1201 Oh 
values stated above and arrive at the percentage increase and hence 
the correction factor for the emergent FUV flux at the metallic¬ 
ity of the galaxy. For spiral galaxies the metallicity at the very 
outskirts of the disc reduces at the most to 0.4-0.6 Zq accord¬ 
ing to available evidence from observations within our Galaxy 
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dCarraro et al.ll200^ : iGenovali et alj|2014h . This also seems to be 
the case for the THINGS sample galaxjos_l[Mous^as_et_t ^2010h 
and the ‘Bluedisk’ sample galaxies dCarlton et. alJlin prep.!) with 
available metallicity gradient measurements. For this metallicity 
the correction factor is only ~10%. This is small compared to the 
total error on the SFR calibration, and therefore we do not apply 
any correction for the reduction in ISM metallicity in the outskirts 
of the spiral galaxies in our sample. 


4 RESULTS AND DISCUSSION 

4.1 Deriving the K-S relations for the different datasets 

In each logarithmic bin of Egas / Egas,atomic we determine and 
plot: (i) the logarithms of the median (median (Log Esfr)), 5 
and 95 percentile Esfrs taking into consideration all regions in¬ 
cluding those with zero or negative fluxes, (ii) the logarithm of 
the average Esfr calculated by summing over the fluxes from all 
regions inc luding those with zero or negative flux es: Log (aver¬ 
age Esfr). Ida Silva. Fumagalli & Krumholzl ( |2014|) show that due 
to the reasons discussed in Section [T] if the SFR being traced is 
low, using standard multiplicative calibrations to convert a mea¬ 
sured FUV flux to SFR will give erroneous results. Even though 
the median flux characterizes the distribution of measured fluxes 
in a robust manner, if the median flux originates from a low SFR 
region the median (Log Esfr) determination is likely to be bi¬ 
ased to a lower value than the true underlying SFR. In the HI 
dominated regime under investigation, most fluxes from individ¬ 
ual regions are in this uncertain regime. A comprehensive ap¬ 
proach of quantifying the relation between surface densities of gas 
and star formation in this regime will involve correcting for the 
bias in the measured SFR using the output f r om si mulations like 
the one by Ida Silva. Fumagalli & Krumholzl l l2014l) . Their simu¬ 
lations show that although the bias in measured SFR is generic, 
the exact amount of the bias depends critically on quantities like 
the total fraction of stars forming in clusters. Such quantities are 
as yet not well constrained from observations for dwarf galax¬ 
ies and in the outskirts of spirals, and therefore it is beyond the 
scope of this paper to determine the corrected statistical distri¬ 
bution of Esfr in the low SFR regime. In this study we there¬ 
fore use the alternative determination Log (average Esfr), for 
which the summed total SFR in the bin always remains above 
log(SFR/(MQ yr~^ )) ~ — 3, where FUV should be a high ly re- 
liable SFR estimator ([^ Silva. Fumagalli & Krumhol^l2014l) . We 
note that we have also checked that the distribution of negative flux 
values is well approximated by a one-sided Gaussian function, and 
that there are no significant tails at large negative values that com¬ 
promise the summation. 

We carry out similar investigations for our different datasets 
covering different galaxy samples and scales. For each dataset 
merely the Log (average Esfr) values are considered for deter¬ 
mining Kennicutt-Schmidt relations. In each case the relation is 
determined using a linear regression fit to the Log (average Esfr) 
and binned Log (Egas,atomic) values, only considering bins with 
a statistically significant number of individual regions in them (~ 
100 or more). It is to be noted that a relation thus determined is 
representative of the ‘average’ trend, and the error on the fit is only 
indicative of the variation in any such ‘average’ relation and is not 
indicative of the scatter that might be present in the Esfrs from 
low Egas,atomic regions. Also the Kennicutt-Schmidt relations de¬ 
termined here are subject to change depending on correction of low 
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Figure 1. The logarithms of the surface densities of SFR and gas 
(atomic+molecular) plotted against each other for 1 kpc sized regions in 
THINGS sample galaxies. The open dark blue squares (in the upper panel) 
represent the median (Log Esfr) (see text) in each 0.3 dex Egas bin and 
the open light blue squares represent the 95 and 5 percentile points (shown 
at y=—8 if negative), with eiTorbars determined by bootstrapping the data 
in each Egas bin. In the upper panel filled red circles represent the Log (av¬ 
erage Esfr) (see text) with associated errorbar. In the upper panel the bold 
line at the h i gh Ega s end is the Eh 2 based K-S law for nearby spirals from 
iLerov et aU i2013ah . Dotted beige lines in the background indicate various 
constant gas depletion timescales. The lower panel shows the fraction of 
points lying below the plotted range including those with negative FUV 
flux in each Egas bin. 

Esfr data for bias as discussed above, and statistically robust mod¬ 
elling to fit the corrected dataset after taking into consideration the 
full scatter in measured Esfrs. A distinctive feature of our study is 
that we select and study the Kennicutt-Schmidt relation for the en¬ 
tire Hl-dominated ISM of the galaxy (discussed below), as opposed 
to concentrating on regions in the disk within some fixed annuli. 

1 kpc scale 

We first show results for THINGS spiral galaxies over 1 kpc sized 
regions in Figure [T] For this dataset we have measurements of 
molecular gas, and can do a direct comparison of our results with 
previous results which used the same sample. At the high Egas end 
where molecular hydrogen dominates and the SFRs in individual 
regions are high, our determinations of both median (Log Esfr) 
and Log (average Esf r) are consistent w ith each other and our 
results match those of jLerov et al.ll20133) very well. At low val¬ 
ued of Egas the scatter in values measured in individual regions 
increases, and the number of regions with very low or negative 
FUV flux and below detection limit increases. We find that the 
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Figure 2. The logarithms of the surface densities of SFR and atomic gas 
plotted against each other for 1 kpc sized regions in THINGS sample galax¬ 
ies. The bold black line represents the best fit Kennicutt-Schmidt relation 
(see text). Other symbols are similar to those in Figure[T] 


K-S r elation for the outer discs of spirals derived by iBigiel et'd] 
(l2010h matches our determination of the median (Log Esfr). For 
low Egas bins, however, the median (Log Esfr) diverges from the 
Log {average Esfr) (middle panel of Figure [TJ. The difference 
between the Log {average Esfr) and median (Log Esfr) reaches 
almost an order-of-magnitude at the lowest gas surface densities. 
The Log {average Esfr) does not exhibit as steep a decline as the 
median (Log Esfr), and the average gas consumption timescale 
remains around few tens of Gyrs at the lowest gas column densities 
we can trace (1 Mq pc“^). 

The results for the FlI dominated regions at 1 kpc scales in 
TFIINGS spirals are shown in Figure and the best fit values to 
the binned Log {average Esfr) - Log (Egas,atomic) values is listed 
in Table |4] We have repeated the study by smoothing all data to 
exactly 1 kpc resolution, and found differences in the results to be 
insignificant. 


400 pc scale 

The results at 400 pc scales in TFIINGS spirals and FIGGS dwarfs 
are shown in Figure 14.11 and the best fit values tabulated in Ta- 
ble|4l An interesting point to note is that for the lowest Egas,atomic 
bin in the THINGS sample plot, a large fraction of regions with 
low or negative Esfrs results in the median (Log Esfr) to de¬ 
crease sharply compared to the value in the immediately higher 
Egas,atomic bin, whereas the decrease in Log {average Esfr) 
compared to the value in the higher Egas,atomic bin is moderate 
and follows the determined Kennicutt-Schmidt relation. As can be 
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Figure 3. The logarithms of the surface densities of SFR and atomic gas 
plotted against each other for 400 pc sized regions in THINGS (upper half) 
and FIGGS (lower half) sample galaxies. The bold black lines represent the 
best fit Kennicutt-Schmidt relations (see text). The grey dashed line in the 
FIGGS sample panel is the best fit relation for disk-averaged values from 
iRovchowdhurv et aljj2014l) . Other symbols are similar to those in Figure^ 
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Table 4. Comparing derived Kennicutt-Schmidt relations 


Dataset 

slope 

intercept 

FIGGS 400pc 

1.5 ±0.1 

-4.36±0.06 

THINGS 400 pc 

1.43±0.07 

-4.35±0.05 

THINGS I kpc 

1.65±0.04 

-4.74±0.03 

Bluedisk control 

1.52±0.06 

-4.18±0.02 

Bluedisk main 

1.71 ±0.08 

-4.64±0.04 


seen from the figure, we restrict the fit to Egas.atomic bins with 
small fraction of regions with noisy Esfr- The apparent mismatch 
between the disc-averaged K-S relation seen in FIGGS galaxies 
dRovchowdhurv et al.ll20l4) and the steeper relation we find here 
at 400 pc scales is explainable by the fact that for the disc-averaged 
study we bias ourselves to ~kpc sized regions which are forming 
stars (the optical disc). The extended HI with little or no star for- 
mation present in most of th e FIGGS galaxies was excluded in the 
iRovchowdhurv et aP ( l2014h study. 

10 kpc scale 

The results for the ‘Bluedisk’ sample spiral galaxies are shown in 
Figure 14.11 and the best fit values tabulated in Table |4] The dif¬ 
ferences in the K-S relation seen in ‘Bluedisk’ main and control 
galaxies is worth noticing. For the main sample the large amounts 
of HI outside the optical disc which have negligible or no star for¬ 
mation going on in them appears to steepen the relation at low gas 
surface density. 

In order to compare the results from the ‘Bluedisk’ spirals 
with that obtained for the nearby THINGS galaxies, we tried 
smoothing the data for THINGS galaxies to 10 kpc scales. But this 
did not result in a trustworthy result, as the FUV maps had to be 
smoothed 20 times or more which resulted in the final FUV maps 
to be dominated by noise except in the very central regions of the 
sample galaxies. We are also aware of the fact that if the THINGS 
galaxies were to be imaged at 10 kpc resolutions, some low level 
highly extended emission may be present which is not picked up at 
their present ~ kpc resolution. Smoothing the data does not alle¬ 
viate this issue. Notwithstanding the existence of these sources of 
uncertainty, a linear fit to the final binned data for THINGS galax¬ 
ies smoothed to 10 kpc resolution matched the values derived for 
the ‘Bluedisk’ control sample well. 

Quantifying the scatter 

In Table [5] we show how the offset between measured Esfrs 
and their scatter in the Esfr - Egas,atomic plane varies with 
Egas,atomic and between datasets. Negative Esfr values are ar¬ 
rived at using the SFR calibration on negative fluxes and have no 
physical meaning, but are used here to be indicative of the scat¬ 
ter. Column (1) gives the centre of the Egas.atomic bin in log, (2) 
the difference in dex between the Log (average Esfr) and me¬ 
dian (Log Esfr) values. For ‘Bluedisk’ galaxies where fluxes are 
averaged over ~10 kpc sized regions (equivalent to summing over 
fluxes from hundreds of 400 pc or 1 kpc sized regions), the differ¬ 
ence is the least and match well considering the errors on the mea¬ 
sured Log (average Esfr)s. For datasets with 400 pc and 1 kpc 
sized regions the difference increases with decreasing Egas,atomic, 
and the offset is more at the same Egas,atomic bin for smaller sized 
regions. Columns (3) through (6) list the 95, 84, 16, and 05 th 
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Figure 4. The logarithms of the surface densities of SFR and atomic gas 
plotted against each other for HI beam sized regions in ‘Bluedisk’ control 
(upper half) and main (lower half) sample galaxies. The bold black lines 
represent the best fit Kennicutt-Schmidt relations (see text). Other symbols 
are similar to those in Figure[T] 
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Table 5. Statistics for Ssfr - atomic plane 


Log Sgas,atomic 

diff. 

percentile SsFR 

in Mq yr 


(Mq pc-2) 

(dex) 

95 

84 

16 

05 

FIGGS 400 pc 

-0.45 

0.71 

6e-05 

8e-06 

-8e-07 

-le-05 

-0.15 

0.66 

7e-05 

9e-06 

-4e-07 

-3e-06 

0.15 

0.62 

3e-04 

8e-05 

2e-06 

-le-06 

0.45 

0.82 

6e-04 

2e-04 

3e-06 

-2e-06 

0.75 

0.59 

2e-03 

6e-04 

3e-05 

9e-06 

1.05 

0.27 

6e-03 

2e-03 

3e-04 

2e-04 

THINGS 400 pc 

-0.15 

undef. 

9e-05 

2e-05 

0 

-le-05 

0.15 

0.60 

4e-04 

9e-05 

0 

-8e-06 

0.45 

0.49 

le-03 

4e-04 

2e-05 

6e-06 

0.75 

0.18 

3e-03 

le-03 

le-04 

4e-05 

1.05 

0.12 

6e-03 

3e-03 

4e-04 

3e-04 

THINGS 1 kpc 

-0.15 

0.57 

4e-05 

8e-06 

-le-06 

—6e-06 

0.15 

0.52 

8e-05 

3e-05 

le-07 

-7e-06 

0.45 

0.53 

4e-04 

le-04 

5e-06 

-8e-06 

0.75 

0.41 

le-03 

6e-04 

4e-05 

8e-06 

1.05 

0.25 

3e-03 

le-03 

2e-04 

7e-05 

1.35 

0.20 

9e-03 

5e-03 

8e-04 

5e-04 

Bluedisk control 

-0.45 

0.18 

5e-05 

3e-05 

le-06 

-2e-05 

-0.15 

0.27 

le-04 

7e-05 

7e-06 

-8e-07 

0.15 

0.20 

4e-04 

2e-04 

2e-05 

9e-06 

0.45 

0.13 

8e-04 

5e-04 

9e-05 

4e-05 

Bluedisk main 

-0.45 

0.17 

2e-05 

le-05 

-2e-06 

-le-05 

-0.15 

0.20 

5e-05 

2e-05 

-3e-06 

-le-05 

0.15 

0.22 

le-04 

6e-05 

—4e-06 

-9e-06 

0.45 

0.19 

3e-04 

2e-04 

3e-05 

le-05 

0.75 

0.20 

2e-03 

7e-04 

2e-04 

le-04 


percentile Esfr value measured in different regions within a re¬ 
spective Egas,atomic bin. The scatter in measured Esfrs obviously 
increases with decreasing Egas,atomic, and decreasing region size. 
But there is no clear indication that the scatter increases for dwarf 
galaxies compared to spirals for the same region size. 

4.2 A universal K-S relation in HI dominated ISM 

In Figure|5]our fits to the Kennicutt-Schmidt relation in the HI dom¬ 
inated ISM of various types of galaxies and at different resolutions, 
are compared to previous determinations of the relation between 
Esfr and Egas in spir als. The linear fit to the relation seen for 
outer discs of spirals bv iBigiel et alj ( l2Q10h is also provided as ref¬ 
erence, although in their study only regions within an annulus with 
inner radius of R25 and outer radius of 2 times R25 were consid¬ 
ered. From the values listed in Table [T] we can see that a K-S type 
relation in HI dominated gas with a power law slope of ~1.5 is a 
good fit to the data for both ^ T* spirals and dwarf galaxies. The 
derived K-S relation is always much more inefficient than that be¬ 
tween molecular/total gas and SFR within the optical disks of L* 
spirals. The amplitude of the relations for different samples agree 
within a factor of ~ 2. 

We have derived the Kennicutt-Schmidt relation in the HI 
dominated ISM of three very different environments - dwarf galax¬ 
ies, outer disks of normal galaxies and outer disks of ‘Bluedisk’ 



Log (^gas ! (1^0 po'^)) 


Figure 5. The resolved Kennicutt-Schmidt relation in the HI dominated 
regime, in relation to other determinations . The black line rep resents the 
Eh 2 based K-S law for nearby spirals from lLerov et alJ42013ah . the brown 
contour indicates the relat ion as seen on 750 pc scales within the optical 
discs of nearby spirals bv lBigiel~ in boosh (the red shaded area in the 
middle right panel of their Figure 8), and the yellow open circles represent 
the median Ssfr in similar Egas bin s in the very outer discs of spirals 
from the THINGS survey measured bv iBigiel et alj iioiot) at 1 kpc scales 
with the scatter indicated by the errorbars. The K-S type relations for the 
different datasets studied in these paper are marked by different colours and 
shown in the range of Egas,atomic within which the respective linear fits 
were done. The 1 a enor on the fits and their overlap are represented by 
the grey shaded area. Dotted beige lines in the background indicate various 
constant gas depletion timescales. 

galaxies - and found them to remarkably similar. One of the prop¬ 
erties which characterizes the diversity of the ISM in the different 
datasets studied, is their metallicity. In Figure]^ we check for any 
variation of HI depletion time with metallicity. For clarity we show 
the results for two datasets mapped at the same resolution: FIGGS 
galaxies and THINGS galaxies mapped at 400 pc resolution. We 
do not have comprehensive measurements of the metallicities for 
all the regions in our sample galaxies, but deduce the expected 
ranges in metallicity sampled by the different datasets using ex¬ 
isting observations and trends. For the FIGGS sample galaxies the 
metallicities are assumed to span the range of values tabulated in 
Table mindful of the fact th at dwarf galaxies do not show ap¬ 
preciable metallicity gradients l lWestmoauette et ^l2013h . For the 
THINGS sample galaxies the lower limit in metallicity is defined 
by the observed flattening in the metallicity gradient b eyond R25 
in nearb y spirals llBe rg et alJl 2013l : lBresolin et al.ll2012h and in our 
Galaxy dGenovali et al.ll2014h . The upper range in metallicity for 
THINGS galaxies is determined by the average measured metal¬ 
licity at the maximum radius from the centre of the sample galax¬ 
ies where CO emission in observ ed, for the galaxies with existing 
metallicity gradient measurements iMoustakas et al.l(l2010ll . We use 
the average of the metallicities (in solar metallicity units) derived 
using the two differ ent metallicity gradient s (following KK04 and 
PT05 methods, see IMoustakas et ani2010L for details). There are 
some Hl-dominated regions at radii smaller than the maximum ex- 
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Figure 6. The vaiiation of HI depletion time with metallicity when mea¬ 
sured at 400 pc resolution. The red point and errorbars represent the value 
for the FIGGS sample. The blue filled points and bold en'orbai's repre¬ 
sent the values for THINGS galaxies (400 pc resolution) in two different 
metallicity ranges. The average value for all Hl-dominated regions in the 
THINGS 400 pc dataset is represented by the blue open square with dashed 
errorbars. 


tent of CO in each sample galaxy disk, but these regions are inter¬ 
arm regions. The metallicity in inter-arm regions of spiral galaxies 
are lower tha n the metallicity me asured from HII regions in the spi¬ 
ral arms (e.g ICedres et al.ll 2012 l) . and we therefore assume that the 
metallicity of these interarm regions are at the most solar. The HI- 
dominated regions from the THINGS dataset are divided into two 
roughly equal rang es in metallicity using t he metallicity gradient 
measurements from iMoustakas et al.l (l2010l) . 

We determine the average HI depletion times by summing up 
the HI column densities and the SFRs from all regions under con¬ 
sideration. The errorbars shown in Figure [^incorporate all uncer¬ 
tainties related to HI, UV and IR flux measurements as well as the 
uncertainty in SFR calibration. For the THINGS galaxies we see 
that when radial distance from the galaxy centre is used as a proxy 
for metallicity of the ISM, the HI depletion time increases with de¬ 
creasing metallicity within spiral galaxies. But we also derive the 
average depletion times for all Hl-dominated regions of THINGS 
galaxies, which is a more appropriate measure to compare with the 
HI depletion time derived for FIGGS galaxies (which do not have 
metallicity gradient measurements). And on comparing the average 
values we notice that the HI depletion times are remarkably uni¬ 
form across the large range in metallicity sampled by dwarf irregu¬ 
lars and spirals. We have checked that THINGS regions mapped at 
1 kpc resolution, and ‘Bluedisk’ control galaxies give comparable 
values for the measured average HI depletion time. 

The amount of dust in the ISM scales with its measured 
metallicity. (i) The rate of H 2 formation and (ii) the shielding 
of UV radiation enabling the physical conditions for H 2 forma¬ 
tion, are both proportional to the amount of dust in the ISM, 
and are therefore expec ted to be proportional to the metallic- 
ity of the ambient ISM ( iKrumholz. McKee & TumlinsonI 1200^ : 


IStemberg et al.l l2014h . There also exist observational evidence 
of the dependence of H 2 -to-HI co nversion on metallicity (e.g. 
IWeltv et^l2012l : IWong et al.ll2013h . The H 2 thus formed is ex¬ 
pected to be the phase correlated with star formation as in this 
phas e cooling and fragmentation of ga s proceeds very efficiently 
(e.g. IKrumholz. Leroy & McKee! |20111) , and as mentioned previ¬ 
ously o bservations seem to confirm this in the central regions of 
spirals jLerov et alj|20133) . One would thus expect that a higher 
efficiency of conversion of HI to H 2 should reflect in a lower gas 
depletion time / more efficient Kennicutt-Schmidt relation. That is 
why our finding that the gas depletion time / Kennicutt-Schmidt 
relation is independent of the metallicity in HI dominated regions 
of star forming galaxies is surprising. Though it must be noted that 
molecular gas formation might not be determined by metallicity 
alone, and factors like the local pressure and volume densities could 
also play a role. 


4.3 Comparison with theoretical models 


In this section, we compare our derived K-S relations with pre¬ 
dicted trends from two wid ely used models of star formation based 
on the followin g papers: ( i)IOstriker, McKee & Lero ^(l201(]t) here¬ 
after OML, (ii) llGumholz. McKee & TumlinsonI 1 2009d) hereafter 
KMT. The OML model assumes that atomic gas has achieved 
two-phase thermal and quasi-hydrostatic equilibrium. In the KMT 
model, hydrostatic balance sets the requirements for the forma¬ 
tion of the CNM phase and the atomic-to-molecular transition, 
except at very low metahicities where it is set by the two-phase 
ther mal equilibrium a t the maximum allowed CNM tempera¬ 
ture jKrumholzl l2013l) . We consider the predictions of the re¬ 
spective models for the low SFR, low metallicity, Hl-dominated 
regime. Therefore for OML mod el we consider equation (22) from 
lOstriker. McKee & LerovI bOlOh . For the KMT model, we con- 
sider t he combination of equations (18) and (20) from iKrumholzl 
(l2013h . 


When comparing the models with our observations four pa¬ 
rameters used in both the OML and KMT models need to be consid¬ 
ered: metallicity, psd, f-w and the clumping factor. For the FIGGS 
galaxies the metallicity is fixed at 0.1 Zq, whereas for the spirals 
from THINGS and ‘Bluedisks’ surveys, the metallicity is varied be¬ 
tween 0.4 to 1 Z 0 . The parameter p^d in the models which refers to 
the density of stellar - |- dark matter is varied betwe en 0.003 to 0.03 
Mq pc“^ (following iKim. Kim & Ostrikeill201ll) . The parameter 
/m is the models which is related to the fraction of n eutral gas in 
the war m phase is varied between 0.05 and 1. From ILerov et al.l 
(l2013bl) we know that HI clumps much less when observed at high 
resolution compared to CO, and in fact the clumping factor for HI 
remains almost constant. As we are only sensitive to HI emis sion, 
we follow the empirical relation given in lLerov et al.l ( 12013ht) and 
adopt a clumping factor of ~ 1.3 at 400 pc and 1 kpc scales (in 
comparison with what will be observed at 100 pc resolution) and ~ 
1.5 at 10 kpc scales. 

Models with the range in parameters given above are shown 
in Figure [7] We see that the predictions of the OML model 
match our observations better than the predictions of the KMT 
model for galaxies of all types. It is interesting that our data 
agrees w ith the OML model, since recent hydrodynamical sim- 
ulations dKim. Kim & Ostrikeill201 ll : iKim. Ostriker & Kimll2013h 
have shown that turbulence driven by supernova feedback regulates 
the thermal pressure in the outer discs of galaxies, and this thermal 
pressure plays a crucial role in the OML models. This may imply 
that the ISM in outer discs of spirals and dwarfs is stabilized by 
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Figure 7. Comparison of the K-S relations determined by us for FIGGS dwarfs (left panel), THINGS spirals (middle panel) and ‘Bluedisk’ galaxies (right 
panel) with OML (blue shaded area) and KMT (red shaded area) models. See text for details of model parameters used. 1 a en'ors on the fits are shown in grey. 


feedback. Recent hydrodynamical simulations of galaxy formation 
have emphasized the prim acy of feedback in de termining how gas 
converts into stars (e.g. see lHopkins et al.ll20l4) . 


5 SUMMARY 

We study the spatially resolved Kennicutt-Schmidt relation be¬ 
tween Esfr and Egas,atomic in the HI dominated ISM in spiral 
and dwarf irregular galaxies across a wide range of spatial scales 
and gas surface densities. We account for the uncertainties associ¬ 
ated with measuring low SFRs within small regions by averaging 
over SFR in fixed bins of HI surface density. For dwarfs to spirals, 
and across different scales, the K-S relation is found to have a slope 
of ~ 1 . 5 . The efficiency with which the cold neutral gas is being 
converted into starts is an order of magnitude more inefficient than 
within the optical discs of spirals. A surprising conclusion from our 
findings in that the mean Kennicutt-Schmidt relation in HI domi¬ 
nated regions is independent of metallicity. Our results favour a 
model of star formation where thermal and pressure equilibrium 
in the outer ISM regulate the rate at which star formation occurs, 
where the thermal pressure in turn is set by supernova feedback. 
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